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ABSTRACT 

The primary goal of this paper is to provide the evidence that can either prove or 
falsify the hypothesis that dark matter in the Galactic halo can clump into stellar-mass 
compact objects. If such objects existed, they would act as lenses to external sources in 
the Magellanic Clouds, giving rise to an observable effect of microlensing. We present 
the results of our search for such events, based on the data from the second phase of 
the OGLE survey (1996-2000) towards the SMC. The data set we used is comprised of 
2.1 million monitored sources distributed over an area of 2.4 square degrees. We found 
only one microlensing event candidate, however its poor quality light curve limited 
our discussion on the exact distance to the lensing object. 

Given a single event, taking the blending (crowding of stars) into account for the 
detection efficiency simulations, and deriving the -ffijT-corrected number of monitored 
stars, the microlensing optical depth is r = (1.55±1.55) x 10 -7 . This result is consistent 
with the expected SMC self-lensing signal, with no need of introducing dark matter 
microlenses. Rejecting the unconvincing event leads to the upper limit on the fraction 
of dark matter in the form of MACHOs to / < 20 per cent for deflectors' masses 
around 0.4 M Q and / < 11 per cent for masses between 0.003 and 0.2 M Q (95 per cent 
confidence limit). Our result indicates that the Milky Way's dark matter is unlikely 
to be clumpy and form compact objects in the sub-solar-mass range. 

Key words: Cosmology: Dark Matter, Gravitational Lensing, Galaxy: Structure, 
Halo, Galaxies: individual: Small Magellanic Cloud 



1 INTRODUCTION 

The Magellanic Clouds are harbours to millions of stars. 
The light of each of these objects can be magnified if an- 
other massive object is close enough to the line-of-sight con- 
necting the observer and a distant star. Paczyhski (1986) 
first realised that with the advent of CCDs, forthcoming 
massive photometric surveys could effectively test the hy- 
pothesis that dark matter in the Galactic halo can clump 

* Based on observations obtained with the 1.3 m Warsaw tele- 
scope at the Las Campanas Observatory of the Carnegie Institu- 
tion of Washington. 
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and form Massive Compact Halo Objects (MACHOs). These 
objects, if they existed, would act as lenses to more dis- 
tant LMC/SMC stars, within the reach of current observ- 
ing facilities. This brilliant, yet simple idea triggered sev- 
eral microlensing programs to emerge. The first detections 
of the microlensing effect were reported by the MACHO 
(Alcock et al. 1993), OGLE (Udalski et al. 1993), EROS 
(Aubourg et al. 1993), MO A (Yock 1998), Angstrom (Kerins 
2008), POINT- AGAPE (Auriere et al. 2001), and WeCaPP 
(Riffeser et al. 2003) microlensing teams. 

For almost two decades, microlensing as an astrophys- 
ical tool has been very successful in finding objects which 
do not emit any or emit little light. The OGLE group alone 
have discovered over 4000 ordinary microlensing events to 
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date. A list of exotic microlensing events includes detection 
of black-holes (e.g., Mao et al. 2002), planets (e.g., Udalski 
et al. 2005, Gaudi et al. 2008, Dong et al. 2009), binary stars 
(Skowron et al. 2007) and also a variety of effects such as the 
parallax (e.g., Smith et al. 2003, Gould et al. 2009), xallarap 
(Assef et al. 2006), etc. 

However, since the microlensing field has evolved into 
a tool nowadays primarily concentrated on finding either 
the most distant or the smallest known planets, Paczyhski's 
original idea has been somewhat forgotten. The primary mo- 
tivation for this paper is to fully explore the existing OGLE 
data to search for microlensing events towards the Magel- 
lanic Clouds. As of 20f0 we have collected approximately 
13 seasons (4 seasons of OGLE-11 and 9 seasons of OGLE- 
III) of data for both the LMC and SMC. In Wyrzykowski 
et al. (2009) (hereafter Paper I) we presented our first es- 
timate of the microlensing optical depth towards the LMC 
from the OGLE-II data. The detection of two events led 
to the optical depth of tlmc = (0.43 ± 0.33) x 10~ 7 . How- 
ever, the MACHO collaboration derived the optical depth 
of tlmc = (1-0 ± 0.3) x 10~ 7 based on their 10 candidates 
(Alcock et al. 2000, Bennett 2005). If this number is com- 
pared to the optical depth for the Galactic halo entirely 
made of MACHOs, r ha i « 4.7 x 10~ 7 (Bennett 2005), it 
gives the fractional contribution of / = tlmc / Thalo ~ 20 
per cent. On the other hand, the EROS collaboration has 
derived tlmc < 0.36 x 10~ 7 , which translates to / < 8 per 
cent only (Tisserand et al. 2007). The OGLE-II estimate 
of / < 10 per cent from Paper I, favours the EROS solu- 
tion but the two detected events are also consistent with 
the expected LMC self-lensing signal. 

The SMC has received somewhat less attention in terms 
of microlensing studies than the LMC. So far, only the EROS 
data were studied systematically and the optical depth of 
tsmc = (1.7 ± 1.7) x 10~ 7 was derived for one microlensing 
event detected in their Bright Stars Sample (Tisserand et al. 
2007). Another study of 5 years of the EROS data gives 
/ < 25% for objects with masses 10~ 7 to IMq (Afonso 
et al. 2003). On the other hand, the MACHO collaboration 
estimated the optical depth to be (2 — 3) x 10~ 7 based on 
their two events (Alcock et al. 1999). The SMC self-lensing 
estimates are in a range of (0.4-1.8) x 10~ 7 from N-body 
simulations by Graff &: Gardiner (1999) and from analytical 
work of Palanque-Delabrouille et al. (1998). 

In this paper we extend our work from Paper I, on 
search for dark matter compact objects in the Galactic halo, 
to an independent SMC data set collected by OGLE during 
its second phase in years 1996-2000. The paper has the fol- 
lowing structure. First, the empirical optical depth estimator 
is described. Then the observational data used in the analy- 
sis are presented in Section 3. Next, in Section 4 the search 
procedure for events is described and its yield presented. The 
detection efficiency of events and the calculation of the op- 
tical depth is discussed in Section 5-7. The paper concludes 
with a discussion of the results. 



2 EXPERIMENTAL OPTICAL DEPTH 

A review of the microlensing-related quantities is given in 
Paczyhski (1986, 1996) and (Gould 2000). In short, the time- 
scale (Einstein radius crossing time) is the only physi- 



Table 2. Error correction coefficients for OGLE-II SMC fields. 



Field 


7 


e 


SMC.SC1 


1.195372 


0.003250 


SMC.SC2 


1.223286 


0.001327 


SMC.SC3 


1.277762 


0.002759 


SMC.SC4 


1.153056 


0.002099 


SMC.SC5 


1.218417 


0.002249 


SMC.SC6 


1.214665 


0.001926 


SMC.SC7 


1.233638 


0.002660 


SMC.SC8 


1.190629 


0.002621 


SMC.SC9 


1.182015 


0.002673 


SMC.SC10 


1.167742 


0.002901 


SMC.SC11 


1.222095 


0.002335 



cal parameter of an event that, in the simplest case, can 
be derived while fitting an observed light curve with the 
microlensing model. The basic point-lens point-source mi- 
crolensing light curve (Paczyhski 1996) is described by 



mag(t) = mag - 2.5 log [f s A(t) + (1 - / s )] ; 



(1) 



where mag(t) is the observed magnitude at a given moment 
of time t, mag is the baseline magnitude (away from the 
peak) and fs is the ratio of the lensed source flux to the 
total flux of stars within the seeing disk (blending param- 
eter). The amplification A(t) then depends on the lensing 
geometry, which changes with time, and can be evaluated 
with the equations 



A(t) 



u(ty + 2 

tt(t) v /u(t) 2 +4 



and u(t) 



«8 + ^, (2) 



where «o is the impact parameter and to is the time of the 
maximum of the peak. 

We calculate the empirical optical depth following 
(Paczyhski 1996) as 



2A»T bs 



iVcv 

E 



(3) 



where T Q b s is the time-span of all observations, N* is the 
total number of monitored stars, iV ev is the total number 
of events, £e, is the time-scale of individual events detected 
with the efficiency of e(<Ei). The non-trivial parts of the 
estimator are the total number of monitored stars and the 
detection efficiency, and they are described in detail in Sec- 
tion 6. 



3 OBSERVATIONAL DATA 

The data used in this work were gathered during the second 
phase of the OGLE survey, which started in 1996 and lasted 
until 2000. The project used its own dedicated Warsaw Tele- 
scope located in the Las Carnpanas Observatory, Chile. The 
details on the instrumentation for the OGLE-II survey can 
be found in Udalski et al. (1997). 

A single 2kx2k CCD chip (pixel size of 0.417 arcsec) 
was operated in the drift-scan mode, giving an actual size for 
each observed held of 2kx8k (14x56 arcmin). There were 11 
fields observed towards the SMC covering in total 2.4 square 
degrees. Their locations are shown in Fig. 1 marked over an 
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Figure 1. Positions of the OGLE-II SMC fields (red rectangles with labels). Also shown are all OGLE— III fields (green squares). The 
two small filled squares show the positions of the HST fields used for our blending determination. Background image credit: ASAS all 
sky survey. 



Table 1. Summary of the OGLE-II SMC fields. 



Field 


RAj2000 


Decj2ooo 


No of stars 


Stellar density 


Density level 








template 


estimated real 


[stars/sq.arcmin] 




SMC-SCI 


0:37:50.89 


-73:29:42.0 


96,220 


158,894 


123 


sparse 


SMC.SC2 


0:40:53.11 


-73:17:29.0 


131,911 


217,668 


168 


sparse 


SMC.SC3 


0:43:57.88 


-73:12:29.0 


196,073 


324,564 


250 


sparse 


SMC.SC4 


0:46:58.59 


-73:07:29.8 


246,137 


406,992 


314 


sparse 


SMC.SC5 


0:50:00.59 


-73:08:46.0 


306,784 


591,366 


391 


dense 


SMC.SC6 


0:53:00.77 


-72:58:40.7 


303,571 


585,153 


387 


dense 


SMC.SC7 


0:55:59.64 


-72:53:32.5 


235,195 


389,400 


300 


sparse 


SMC.SC8 


0:58:57.59 


-72:39:30.7 


188,159 


311,264 


240 


sparse 


SMC.SC9 


1:01:54.61 


-72:32:32.5 


159,055 


263,061 


203 


sparse 


SMC.SC10 


1:04:50.47 


-72:24:47.2 


132,781 


219,236 


169 


sparse 


SMC.SCI 1 


1:07:45.40 


-72:39:32.3 


110,230 


182,010 


141 


sparse 


total 






2,106,303 


3,649,608 







Note: Coordinates point to the centre of the field, each being 14' X 56'. Number of "good" objects in the template is provided (TV > 80 
and (I) < 21.0 mag) together with the estimated number of real monitored stars (see Section 6). Stellar density in number of stars per 
square arc minute was used to classify fields into dense or sparse classes with the threshold of 300 stars/sq.arcmin. 
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Table 3. Selection criteria for search for microlcnsing events in the OGLE-II data 



Cut no. No. of objects left 

Selection of "good" objects N > 80, (J) < 21.0 mag 2,106,303 

1 Significant bump over baseline ^ ' <r; > 30.0 3,097 

peak 

2 "Bumper" cut* (I) > 19.0 mag, (V - I) > 0.5 mag 2,295 

x 2 — x 2 

3 Microlensing fit better than constant line fit 2 — > 140 335 



4 Number of points at the peak* Npeak > 6 311 



5 Microlensing fit better than supernova fit Xsn m * n (x 1X^4) 



201 



6 Peak within the data span 466 sg to «S 1874 192 
[HJD-2450000] 

7 Blended fit converged < / s < 1.2 44 

8 Conditions on goodness of microlensing fit — sg 2.3 and x ^<p' ak — ^ 5 2(1)* 
(global and at the peak) 

9 Time-scale cut 1 < *e < 500 2(1)* 

M 

10 Impact parameter cut < uq ^ 1 2(1)* 



t magnitudes as in the field SMC-SC1 (shifted according to the position of the centre of Red Clump) 
tin the range of to ± Kb 

* there was in fact only one event occurring, whose flux was detected on two template objects (see text) 



Table 4. OGLE-II database stars in the SMC on which the flux from the microlensing candidate was detected. 



field 


DB 


image 


baseline 




star id 


x [pix] 


y [pix] 


I [mag] 


V [mag] 


SMC.SC7 


193726 


1517.87 


6448.29 


15.34±0.01 


16.48±0.02 


SMC.SC7 


351179 


1515.84 


6448.42 


19.7±0.2 


21.9±0.5 



image of the SMC obtained by the ASAS survey(Pojmanski 
1997). Table 1 lists all the fields with the coordinates of 
their centres, number of good objects in /—band, blending- 
corrected number of stars (see Section 6) and blending den- 
sity group. By "good" we mean all objects having at least 
80 observations during the entire time span of the OGLE- 
II (from about 22 to 27 per cent of all collected frames of 
a field) and mean magnitude brighter than 21.0 mag. The 
limiting magnitude was chosen at the peak of the observed 
luminosity functions. Similarly as in Paper I, fields were di- 
vided into dense and sparse according the density of objects. 
Here we set the density boundary at 300 "good" stars per 
sq. arcmin, slightly less than in Paper I, in order to diversify 
the density levels in the SMC. The density levels correspond 
to respective HST images which were used for deriving the 
blending distributions (see Section 6). 

The OGLE-II observations of the SMC began in Jan- 



uary 1997 (HJD=2450466) with the monitoring of the SC5 
and SC6 fields only. Nearly half a year later (HJD=2450621) 
the remaining nine fields were added to the observing queue. 
All eleven fields were monitored continuously until Novem- 
ber 2000 (HJD=2451874) yielding about 300 frames per field 
in /-band. Additionally, there were 30 to 45 frames per field 
collected in V-band. On average each field was observed ev- 
ery third night in /-band and every 8-th night in 1^-band, 
with the mean seeing of 1.36 arcsec in / and 1.39 arcsec in 
V. Raw images were de-biased and flat-field corrected "on- 
the-fly" . The photometric pipeline was based on the Differ- 
ence Image Analysis method (DIA, Wozniak 2000; Alard & 
Lupton 1998). The template images for DIA were created by 
stacking the best quality images, resulting in images with the 
seeing of about 1.1 arcsec. Photometric databases were de- 
signed for both pass-bands as described in Szymahski (2005) 
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Figure 2. Astrometric analysis of the OGLE-SMC-01 microlens- 
ing event is shown. Top panel shows the light curve of the event, 
divided into parts: rising (green), maximum (red) and falling 
(blue). Bottom panel shows the 8x8 pixel area from the OGLE-II 
template, centred close to the event position. The range of grey 
levels starts from the lack of flux (white) and increase towards 
dark levels with the increasing flux. The centroids of the bright 
and faint stars from the OGLE-II catalogue, on which the event 
was detected, are marked with the plus and star, respectively. 
With the colour coding from the top panel, we show the posi- 
tions of the magnified light measured on the subtracted images. 
The X symbol shows the mean position derived from the indi- 
vidual positions. This event location is inconsistent with any of 
the two detected stars, meaning that there must be yet another, 
third object that was microlensed. 

and are available on-line 1 . The photometry of the template 
images contain about 2.1 million objects in I— and V— band 
colours suitable for our study. The search for microlensing 
events and determination of the detection efficiency were 
performed using the I— band data only, as these were far 
more numerous and were sampled more frequently as com- 
pared to the V— band light curves. In this paper we also 
occasionally use additional data collected during the third 
phase of OGLE (OGLE-III), years 2001-2009. The detailed 
description of the OGLE-III instrumentation and photomet- 
ric pipeline can be found in Udalski (2003). The search for 
microlensing events in the OGLE-III SMC data will be pre- 
sented separately in a forthcoming paper (Wyrzykowski et 
al. in preparation). 

Before we start the search procedure, all /-band pho- 
tometric measurements have to have their error-bars cor- 
rected, since DIA is known to underestimate the photomet- 
ric error-bars (Wozniak 2000). The correction technique is 
detailed in Paper I. In short, the method compares the in- 
trinsic error-weighted rms of constant stars in each field with 
their mean error. The error-bars are corrected with the for- 
mula 



1 http://ogledb.astrouw.edu.pl 



AI cor = ^( 7 A/)2 + e 2, (4) 

where AI cor is the corrected error-bar, AI is the original 
error-bar returned by the photometry pipeline, and 7 and 
e are the correction coefficients. They are derived for each 
SMC field and are presented in Table 2. The mean 7 and e 
for all SMC fields were 1.20715 and 0.002436, respectively. 



4 SEARCH PROCEDURE 

The main objective of this paper is to find rare microlensing 
events amongst millions of SMC stars. The task requires a 
number of steps to iteratively remove all unwanted, contam- 
inating light curves. 

Our search procedure begins by pulling out light curves 
from the database with more than 80 epochs and the mean I- 
band magnitude brighter than 21.0 mag (Cut 0). This sam- 
ple of 2.1 million light curves becomes our All Stars Sample. 
The following steps along with the number of objects left af- 
ter each cut are detailed in Table 3. We use the same search 
pipeline as in Paper I, however, the exact parameters of the 
cuts were fine-tuned independently for OGLE-II SMC data 
with the Monte Carlo simulations. One of the main differ- 
ences is in Cut 0, where the depth of the search has changed 
from 20.4 mag (LMC) to 21.0 mag (SMC) due to a lessen 
crowding and greater distance modulus. Another important 
change, as compared to Paper I, is the definition of the "blue 
bumper" region (Cut 2). The magnitude limit is chosen to 
be 0.5 mag fainter, than for the LMC, due to larger distance 
modulus of the SMC. 

A potential source of contamination in a sample of mi- 
crolensing events can be blue bumpers and supernovae. Both 
these classes are characterised by similar light curves to that 
of microlensing, especially if a light curve is of low photo- 
metric quality or sparsely sampled. Blue bumpers can be ef- 
fectively removed by using Cut 2, as they occupy the bright 
blue end of the main sequence. Also blue bumpers are known 
to exhibit several bumps on a time-scale of a few years. This 
is why we make use of the OGLE-III data, giving us a full 
span of ~ 13 years, to understand the nature of each mi- 
crolensing event candidate. The contamination with super- 
novae in the OGLE-II is very limited due to the fact that 
we observed a very small central region of the SMC. There 
should be statistically 1 (4) SNe peaking above / < 20 (21) 
mag in our data. In Cut 5, we compare the goodness of fit 
of the microlensing and supernova models, removing ~ 100 
light curves with an asymmetric bump. 



5 SEARCH RESULTS 

Our search pipeline has returned two microlensing event 
candidates. The close inspection of their positions, how- 
ever, revealed they are separated from each other by 2 pixels 
(~ 0.9 arcsec) only. Moreover, the moment of the maximum 
brightness was exactly the same in both light curves. Table 
4 presents detailed information regarding both objects, in- 
cluding their positions and magnitudes. The probability for 
two independent microlensing events occurring at the same 
moment of time in such proximity is close to zero. We con- 
clude, therefore, that both these events are the image of the 
same single microlensing event into two nearby stars. 
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To pinpoint the unbiased position of the event we em- 
ployed DIA. Since the template image is matched astromet- 
rically and photometrically to each image in a series, what 
is left on subtracted images are only variable objects (with 
either positive or negative fluxes) and systematic problems. 
All constant objects are subtracted out. Fig. 2 shows as- 
trometric measurements of the centroid of the differential 
flux at various stages of the event (indicated with different 
colours). It is evident that the most amplified data points 
from the peak of the event (also with the most accurate as- 
trometry) are concentrated away from the centroids of both 
"ghost" microlensing events (marked with a plus and a star). 
This indicates that the microlensed flux is not linked with 
neither of these objects and there must be yet another star 
hidden in the wings of the brighter object. This object is 
not present on the template. We repeated the DIA analy- 
sis with a prior knowledge of the exact source location on 
the CCD chip (x,y = 1516.74,6447.66) and measured the 
flux at that position. This microlensing event was dubbed 
OGLE-SMC-01. 



5.1 OGLE-SMC-01 

The OGLE-SMC-01 microlensing event (RA, Dec) = 
(0:56:45.89, -72:37:19.8) appeared on a bright blend of 
at least two stars with the total I = 15.340 ± 0.003 and 
(V - 1) = 0.759 ± 0.009 in the SMC_SC7 field. The J-band 
light curve of the event, measured at the difference flux po- 
sition, is shown in Fig. 3. 
MA CHO data of the event. 

To interpret the nature of this event we also used the SMC 
data collected by the MACHO group, which are publically 
available 2 . In the MACHO database we found a single ob- 
ject (ID 207.16370.17) with the photometry available in 
MACHO-B and MACHO-B bands. Because of somewhat 
larger seeing disk in the MACHO data, as compared to 
OGLE, only one object was visible in the MACHO images, 
corresponding to a blend including both sources detected 
by OGLE. The original MACHO light curves showed a high 
level of noise and variability over entire light curve in both 
bands. To increase the signal-to-noise ratio in these light 
curves we re-reduced the original MACHO images using DIA 
in a similar manner as in Kozlowski et al. (2007). The red 
images suffered from a number of bad columns in the vicin- 
ity of our object, making it impossible for DIA to produce a 
reasonable light curve. We were able to perform the photom- 
etry in the blue channel (Fig. 3, bottom panel). As expected 
the DIA reduction lowered the amount of noise, as com- 
pared to the original MACHO reduction, clearly revealing a 
periodic variability in the baseline. With the period of 369 
days (close to 1 year) this could be some kind of artefact, 
such as differential refraction. The nearby stars with simi- 
lar colour and brightness, however, showed no variability of 
that kind. On the other hand, the lack of any variability in 
the OGLE-III V-band data and in the EROS B-band data 
(P. Tisserand, private communication) leave us at a loss for 
explanation. In the modelling we treated the variability as if 
it was coming from the blended star, and included it in the 
microlensing fit as described in Wyrzykowski et al. (2006). 



http: / /wwwmacho.mcmaster.ca/Data/MachoData.html 
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Figure 3. Light curves of the OGLE-SMC-01 microlensing event 
candidate. Top: OGLE— II at the peak with the best model fit 
(seven-parameter) and its residuals. Middle: full span of OGLE- 
II and available OGLE-III data for that star showing no other 
"bumps" over 11 years. Bottom: full span MACHO- B data trans- 
formed to the OGLE-V along with the seven-parameter mi- 
crolensing model and its residuals. Period of variations is close 
to 1 year and is not present in OGLE and EROS data. 



In order to obtain the location of the source on a CMD, 
we transformed the instrumental MACHO B-band photom- 
etry to the standard V-band. It was done by comparing the 
photometry of the field stars from OGLE I- and U-band, 
and the MACHO B-band. We found the following colour 
transformation: (B-I) = (1. 14694±0.00277)(V -I) + const, 
where const is an artefact from instrumental magnitudes de- 
rived for MACHO data. Using this formula we were able to 
transform the MACHO B-band light curve to the standard 
V-band. Fig. 3 shows the data for the event in OGLE /-band 
and I/-band. 

Microlensing model fit. 

To find the best set of parameters describing the event, we 
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Table 5. Parameters of the microlensing model fits to the event OGLE-SMC-01. Four- and five-parameter fits were performed on the 
OGLE-II /-band data solely. The seven-parameter fit used also the V-band measurements (obtained from transformed MACHO--B data) 
and included the baseline variability model with period of I yr. 



OGLE-SMC-01 

parameter four-parameter fit five-parameter fit seven-parameter fit 



to 


1385.7 


±1.4 


1385.9 


±1.4 


1384.10 


±0.78 


*E 


31.5 


±1.2 


65.0 


±21.8 


89.7 


±14.0 


«o 


1.960 


±0.022 


0.6593 
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Figure 4. Left: CMD of the SMC region around the OGLE-II candidate event with OGLE (red background) and HST (black dots) 
measurements. The open red circle marks the baseline position of the OGLE-SMC-01 event, whereas the full dot indicates the position 
of the lensed source derived from the microlensing model of the event. Events detected by MACHO are also shown: MACHO-97-SMC- 
1 (blue triangle) and MACHO-98-SMC-1 (green star). Dashed lines delimit the "Blue Bumper" exclusion region used in our search 
pipeline. Right: Positions of OGLE and MACHO candidate events on the OGLE-III SMC Red Clump stars count map with grey levels 
corresponding to counts per square arcminute. Yellow contours follow Red Clump density from 100 to 250 stars per sq.arcmin in intervals 
of 50 stars per sq.arcmin. 



fitted a microlensing model to the available light curves in a 
number of ways (Table 5). A simple microlensing event can 
be successfully described with five parameters (Eqns 2 and 
3). Our first fit was done to the OGLE /-band data solely, 
with the blending parameter /g fixed to 1, meaning that we 
assume that all the light comes from the microlensed source 



only and there are no other sources of light contributing 
to the overall flux. Such a four-parameter fit almost always 
converges and returns good first guess parameters for the 
five-parameter fit. Next, we free the blending parameter and 
fit the light curve to obtain almost identical goodness-of- 
fit values as those returned by the four-parameter fit. The 
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derived blending indicates that only 9 per cent of the light 
comes from the source. Such a drop from 100 to 9 per cent 
for the blending parameter, and basically identical goodness 
of the fit, gives an overall impression of the low quality of 
the light curve. 

To improve the situation we performed a simultaneous 
fit to the OGLE J-band and MACHO B-band data. Since 
microlensing is achromatic, the parameters describing the 
geometry of the event (to, t&, Uq) stay the same for each data 
set. The model has only two additional parameters for an 
additional set of data: the blending parameter and baseline 
magnitude. In this model we also included the variability 
component to accommodate the variations present in the 
MACHO data, approximating their shape with 2 harmonics 
with period of P = 369 ± 2 days and T min = 269.56 ± 0.01 
for the epoch of the minimum. Such a seven-parameter fit 
gives the best constraints for the parameters. The % 2 of this 
fit is the smallest from all three fits. 

From the seven-parameter fit we obtained the colour 
and brightness of the source, and their formal errors as fol- 
lows: (V-I)s = 0.46±0.03 mag and I s = 18.67±0.01 mag. 
Note that the error-bar on the magnitude of the source is 
much smaller than calculated from the values and error- 
bars for the fitted parameters. This can be achieved because 
blending parameters of the fit are highly correlated. The 
microlensing model fit to the transformed V-band data, in- 
cluding the baseline's variability, is also shown in Fig. 3. 
We should emphasise the fact, that the model with variabil- 
ity assumes the shape of the baseline variability does not 
change in time and thus can be predicted for the duration 
of the event. The fit to the baseline (Fig. 3) seems reason- 
able. However, noise in the data can hide some irregulari- 
ties of the periodicity. If the above assumption is not valid, 
namely, that the baseline variability is invariant throughout 
the event, there is no other way of removing the variabil- 
ity, and thus the colour of the source can not be derived 
correctly. 

CMD location of the source. 

Fig. 4 shows the derived location of the source (filled red 
dot) and the blend (open red dot) on the colour-magnitude 
diagram, along with two other SMC microlensing events 
MACHO-97-SMC-1 and MACHO-98-SMC-1. The position 
of OGLE-SMC-01 indicates it lies on the edge of the re- 
gion we have excluded in our microlensing candidate search 
pipeline due to possible contamination of "Blue Bumpers" . 
It means, if it was indeed a microlensing event and was not 
blended with any other star, it would have been rejected by 
our search pipeline. On the other hand, the colour and mag- 
nitude of the source were derived because the microlensing 
fit indicated that severe blending was present in both bands. 
If this event is not due to microlensing, we do not have any 
other means to de-blend the observed fluxes, thus the real 
CMD position of the source star in that scenario remains 
unknown. High-resolution imaging might help in disentan- 
gling the true nature of the source star, however, proximity 
of 15.3 magnitude star could be a serious obstacle. 

The accumulation of problems related with this event 
(blending with a very bright star, artefactual variability in 
the MACHO data, etc.), may cast some doubts on the re- 
liability of the colour and brightness derived for the source 
star. Also, the derived colour relies strongly on the accuracy 
of the model of the variability and its predictability during 



the event. Therefore, with the available information, we are 
unable to conclude firmly whether the event is caused by 
microlensing and where its source star lies. 
Location of the blend. 

Another interesting thing to notice is the location of the 
blend object, as it lies in the not very populated region of 
the CMD. A plausible explanation for this could be that 
the bright star does not belong to the SMC and is located 
in the foreground, e.g., a main-sequence star at a distance 
between 1 and 6 kpc. Another reason for this somewhat 
unusual location could be that apart from the lensed source, 
which is bluer, there is another blue component present in 
the blend and the bright star itself is an SMC red giant. 



5.2 Notes on SMC events discovered by MACHO 
and EROS groups 

During OGLE-II observations of the SMC two other events 
happened, which were detected by the MACHO group and 
also reported by the EROS collaboration. The peak of the 
event MACHO-97-SMC-1 (Alcock et al. 1997) occurred well 
before the OGLE-II observations started, therefore only a 
portion of the fading part of the light curve is present in the 
OGLE-II data (star SMC_SC8.207700). It was not enough 
to be detected as a full event by our search pipeline (mainly 
due to Cut 6 on to, see Table 3). We note in the remain- 
ing OGLE-II data, as well as in the entire 8 years of the 
OGLE-III data for this star, there is no other "outburst" 
visible. Lack of secondary deviations over many years and 
the presence of parallax effects in its light curve (Assef et al. 
2006) confirms this is a genuine microlensing event, despite 
it being located in the "Blue Bumper" region on the CMD 
(see Fig. 4). 

The next event detected a year later, MACHO-98-SMC- 
1, aka EROS2-SMC-1 (Alcock et al. 1999, EROS Collabora- 
tion et al. 1998), was an obvious caustic-crossing binary lens 
event, which was observed by many groups and analysed 
thoroughly (e.g.,Afonso et al. 2000). However, the source 
star of this event was too faint to be detected on the OGLE- 
II template images (Udalski et al. 1997), therefore this event 
is not present in the database studied here. However, with 
a baseline of about 22 mag it would not contribute to our 
optical depth calculations due to Cut on the magnitude. 



6 BLENDING AND DETECTION 
EFFICIENCY 

In this study we applied the same technique as in Paper I to 
deal with blending in the SMC fields. In the Hubble Space 
Telescope (HST) archive 3 we identified two images taken in 
the F814W filter (the closest to OGLE's I band) with ex- 
posure times longer than 10 minutes, located in the dense 
and sparse parts of the SMC. Then, for each of the images 
we cross-matched all stars with visible objects on the cor- 
responding OGLE image. This allowed us to derive mean 
blending distributions for dense and sparse fields and the 
distribution of the number of stars in each OGLE object. 



3 http://archive.stsci.edu/hst/ 
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Figure 5. Distributions of blending fractions for stars from the 
analysis of the archival high-resolution HST images of parts of 
the OGLE SMC fields. 



Figs. 5 and 6 show the distributions derived in three mag- 
nitude bins for dense and sparse OGLE fields (see Table 1). 
When compared with similar distributions for the LMC, it is 
apparent that crowding in the SMC fields is smaller. How- 
ever, the blending is still present at some small level and 
affects the way we see the SMC with OGLE instruments. 

Similarly, as in Paper I, we estimated the real number of 
monitored stars in the SMC fields using the distribution of 
the number of stars in a single OGLE object. The observed 
OGLE luminosity function was convolved with the blending 
distribution and total number of stars was calculated. The 
estimated numbers of monitored stars in each field are pro- 
vided in Table 1. On average, the correction factor for dense 
fields was about 1.9 and for sparse fields about 1.7. The to- 
tal number of monitored stars was estimated to be about 
3.6 million, compared to about 2.1 million objects detected 
on the OGLE-II template images in Cut 0. 

Derived blending distributions were then used in de- 
riving a detection efficiency for microlensing events. Details 
of the simulations and the method are provided in Paper I. 
The efficiency of event detection was obtained for each SMC 
field separately for time scales ranging from 1 to 1000 days. 
Because our search algorithm was not sensitive to binary 
lens events and other exotic events, the efficiencies were ad- 
ditionally corrected by a factor of 0.9 to compensate for an 
estimated fraction of around 10 per cent of all events being 
non-single. 

Apart from the All Stars Sample (with the magnitude 
cut as in the Cut 0.), the efficiencies were also derived for the 
Bright Sample of stars, defined as in Paper I, i.e., 1 mag be- 
low the centre of the Red Clump (usually around 19.5 mag). 
Fig. 7 shows an efficiency curve for the sparse field SMC_SC7 
(the one with OGLE-SMC-01) for All and Bright Stars sam- 
ples, with and without blending taken into account. Efficien- 
cies for the remaining fields varied within around 25 per cent 
from these curves, depending on the density level and the 
time span of observations. As expected, the detection efn- 



Figure 6. Distribution of the number of HST stars in one OGLE- 
II object obtained for dense and sparse fields in the SMC. 

ciency for bright sources is significantly higher than that for 
all sources. 

There is also a strong dependence on taking blending 
into account. If blending is neglected, the efficiency rises by 
approximatively 60 per cent, but the number of monitored 
stars is then the same as the number of objects visible on 
the template. However, the smaller number of monitored 
stars (2.1 million compared to 3.6 million after correcting 
for blending, i.e., 70 per cent drop) does not fully compen- 
sate the rise in the detection efficiency. This agrees with our 
conclusions from Paper I that blending should be scrupu- 
lously taken into account when measuring the microlensing 
optical depth. 

Fig. 8 presents the luminosity functions of dense and 
sparse OGLE-II SMC fields (SC6 and SC9, respectively) 
for which blending was derived using HST archival images. 
Upper panel shows functions as observed by OGLE-II, while 
the lower shows estimated underlying luminosity functions 
recovered using blending distributions from the HST. Blend- 
ing is somewhat larger in the denser field, however, the dif- 
ference here is much smaller than in the LMC (see Paper 
I)- 



7 OPTICAL DEPTH ESTIMATE 

Even though the nature of the only candidate event we de- 
tected is unclear, we attempted to determine the optical 
depth of the single event with respect to the All Stars Sam- 
ple. In the Equation 3 we used T b s = 1408 days for the 
time-span of all observations, N, = 3.65 x 10 6 for the esti- 
mated total number of monitored stars (see Section 6) and 
iV cv = 1 for the total number of events. 

As for the time-scale of the event, we first used the 
value obtained in the fit to the OGLE-II /-band data solely 
(t E = 65.0 ± 21.8) and used the efficiency corresponding to 
that time-scale. Table 6 gathers all the values used in the 
calculations. The error in r was calculated with the formula 
given by Han & Gould (1995). 
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Figure 7. Detection efficiencies for OGLE-II SMC.SC7 field for 
All (thin lines) and Bright (thick lines) Stars Samples, with blend- 
ing included (solid lines) and neglected (dashed lines). 



Table 6. The optical depth for the SMC OGLE-II candidate 
event for time-scales derived for different data sets. 



event 


*E 
[days] 
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efficiency not corrected for binary events 
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0.142281 
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(OGLE I data solely) 
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1.74 ± 1.74 


efficiency corrected for binary events 




OGLE-SMC-01 


65.0 ±21.8 
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1.55 ± 1.55 


(OGLE I data solely) 








OGLE-SMC-01 

(OGLE / + MACHO B) 


89.7±&g 


0.141850 


1.93 ± 1.93 



For the original detection efficiency (i.e., not corrected 
for binary events), the optical depth was derived to be 
TSMC-02 = (1.40±1.40) x 10 -7 . If the efficiency is corrected 
for non-detectability of binary lenses, the optical depth be- 
comes TSMC-02 = (1.55 ± 1.55) x 10~ 7 . 

We also calculated the optical depth for the time-scale 
derived for the combination of the OGLE I and MACHO 
B data, which was slightly larger than for the OGLE data 
only, £e = 89.7l}g'g, but still within the la range of the 
former. Because the events' search pipeline and the detection 
efficiency determination did not use MACHO data, we could 
only estimate efficiency, assuming it was similar to the one 
for the OGLE I data solely. After correcting the efficiency 
for missed binary lenses, the optical depth was tsmc-02 = 
(1.93 ±1.93) x 1CT 7 . 




14 16 18 20 22 24 

magnitude 



Figure 8. Top: observed luminosity functions in two OGLE- 
II SMC fields, where the dotted and solid lines represent sparse 
(SC9) and dense (SC6) fields, respectively. The vertical line shows 
the cut-off at / = 21 mag. Bottom: luminosity functions for the 
same fields recovered after applying the blending correction (see 
Section 6). Their shapes follow the prototype LFs of two HST 
fields (dense and sparse) used for the blending correction (thick 
dash-dotted lines). 

8 DISCUSSION 

8.1 On the nature of the event 

The search for microlensing events in the OGLE-II data 
towards the SMC yielded a discovery of a single, rather weak 
candidate event. The nature of the event is not obvious due 
to numerous complications in the available data (see Section 
5.1). 

Variable star scenario. 

The first scenario is that the event is some kind of a vari- 
able star with an occasional outburst. Astrometry of the 
observed additional flux (Fig. 2) indicates that the source 
of the outburst was severely blended with surrounding stars 
and was located in the wings of the bright star. Therefore 
it was impossible to derive its real magnitude and colour 
and, by placing it on the CMD, to link the event to any of 
known outburst-like variables. Judging just from the sym- 
metric shape of the light curve it could be a "Blue Bumper" 
as these are known to resemble microlensing events (Alcock 
et al. 1997). On the other hand, the bumpers tend to re- 
peat their outbursts after several years (e.g., Tisserand et al. 
2007), which we do not see in the available data. The base- 
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line of the event observed by OGLE-II and OGLE-III re- 
mained constant for over 11 years. The periodic variability 
present in the MACHO data was not observed in any other 
data sets and its period close to 1 year indicates it must be 
some sort of observational artefact. Apart from that, MA- 
CHO data do not show any other additional outbursts, so 
when combined with OGLE data the event remains singular 
for more than 15 years of continuous observations. 
Microlensing scenario. 

The microlensing scenario has also its weaknesses. Because 
of the severe blending with a bright star, the deviation in 
the light curve is rather small (~ 0.08 mag). Nevertheless, 
fitting to both OGLE and MACHO data (seven-parameter 
microlensing model) resulted in the source's brightness and 
colour being derived and the source's location placed on the 
CMD (Fig. 4). 

The fact that the source lies in a sparsely populated 
region of the diagram between the main sequence branch 
and red clump suggests the source could be rather nearby, 
located in the foreground halo or the disk of our Galaxy. 
However, lensing in such case is extremely unlikely, as its 
optical depth is of order of 10 -9 (e.g., Han 2008), i.e., 2 or- 
ders of magnitude smaller than that for self-lensing of the 
SMC. On the other hand, the source could be a binary in 
the SMC, with components of different colours blended and 
close enough for the microlensing not to resolve them, as no 
binary-source features are present in the light-curve. 

Moreover, if the lensed source was located at the far 
end of the SMC it would suffer from high internal extinc- 
tion, thus its position on the CMD would move from the 
main sequence towards the less populated region. The cen- 
tral location of the event on the SMC map coincides with the 
region where SMC is the thickest. However, the internal ex- 
tinction derived by Subramanian & Subramaniam (2009) for 
the region of the OGLE-SMC-01 was E(V - 1) w 0.08 mag, 
far too small to explain the unusual position of the source. 
On the other hand, the authors notice that their estimate for 
the extinction could be underestimated as the extinction in 
the case of the LMC, derived using main sequence OB stars 
as tracers (Harris et al. 1997), differs by about 0.2 mag. 
This could indicate some significant differences in the spa- 
tial distributions of different tracers within the Magellanic 
Clouds. Given the source, if reddened, belongs to the O, B 
or A main sequence stars, the internal extinction could reach 
up to ~ 0.3 mag. That would be consistent with the source 
being located in the far end of the Cloud and self-lensed 
by an object from within the SMC. Self-lensing events were 
already suggested to suffer from higher extinction in (Zhao 
1999). 

Interestingly, both other known SMC events from the 
period of the OGLE-II, namely MACHO-97-SMC-1 and 
MACHO-98-SMC-1, were considered as caused by self- 
lensing (Sahu & Sahu 1998, Afonso et al. 2000,Assef et al. 
2006). 

The self-lensing optical depth of the SMC was estimated 
by Graff & Gardiner (1999) in the N-body simulations of the 
SMC events. They concluded the self-lensing contribution to 
the overall r over the EROS-2 fields was tsl ~ 0.4 x 10~ 7 , 
with a maximum value at the very centre of the SMC of 
t = 1.6 x 10~ 7 . These estimates, however, were based on 
rather small line-of-sight depth of the SMC (about 2kpc) 
and are probably underestimated in view of recent results 



from Subramanian & Subramaniam (2009), who derived the 
depth of the SMC bar to about 5kpc. Moreover, because 
OGLE-II covers only the very central parts of the SMC, it is 
more prone to see self-lensing events than EROS, which cov- 
ered a much larger area. The long time-scale of the OGLE- 
SMC-01 event is also in agreement with the expected time- 
scale for self-lensing events of about 100 days for a mean 
lens mass of 0.35 M©. The time-scale of the event is also 
similar to those of both MACHO events, suspected to be 
due to self-lensing. Therefore, we are inclined to conclude 
that the most likely scenario for the OGLE-SMC-01 event 
is the self-lensing in the SMC. 

8.2 On the optical depth towards the SMC and 
the upper limit on MACHOs in the Galactic 
halo 

Values of the optical depth derived for different time-scales 
of the event, obtained using either OGLE data or a combi- 
nation of OGLE and MACHO data, vary in range between 
1.4 and 1.9 10~ 7 . They agree with estimates for self-lensing 
from Graff & Gardiner (1999) for the central parts of the 
SMC, but their statistical significance is very difficult to 
judge because of just one contributing event. They are, how- 
ever, still in agreement with EROS's calculation for the sole 
event MACHO-97-SMC-1, t S mc = (1.7 ± 1.7) x 10~ 7 (Tis- 
serand et al. 2007). Based on their two events, the MACHO 
collaboration estimated r — (2 — 3) x 10~ 7 (Alcock et al. 
1999), which, again, is in agreement with our result. 

Yet, we must emphasise that the value of the optical 
depth obtained just for one single event is statistically not 
very significant. If converted to a fraction of MACHOs in the 
Galactic halo it would mean that between 25 and 40 per cent 
of halo's mass is contained in dark matter compact objects. 
Nonetheless, that would be valid only when the event we 
found was assumed to be caused by MACHO lensing and the 
modelled optical depth was rf ota} = 5.1 x 10~ 7 for a Galaxy 
halo entirely constructed of MACHOs, according to model S 
of Alcock et al. (2000). Furthermore, if we speculate that all 
three events detected by OGLE-II towards both Magellanic 
Clouds were due to the MACHOs, that would mean that 
MACHOs compose about 15 per cent of the total mass of 
the Galactic halo with mean mass of 0.7 M©. This is shown 
as a box in Fig. 9. 

On the other hand, if all three events detected towards 
the SMC and LMC were caused by the regular luminous stel- 
lar component of each Cloud (self-lensing), there is a nil re- 
sult for candidates for the MACHO lensing among 15 million 
monitored stars from both Clouds. That could be translated 
into an upper limit for MACHO fraction in the halo using 
our detection efficiency and time-scale distribution for halo 
lenses from the model S. We took the mean observing time 
(both Clouds were observed for the nearly the same time) 
and we used mean detection efficiencies to derive the total 
number of expected events in LMC and SMC. The efficiency 
and number of monitored stars in the SMC were adjusted 
to be consistent with the magnitude cut applied in the LMC 
study (20.4 mag). We estimated the number of events ex- 
pected to be observed by OGLE-II in the SMC and LMC 
as caused by MACHO lenses only. Following the suggestion 
of Moniez (2010), we assumed all 3 detected events are the 
expected self-lensing signal (N = 3, B = 3) and used Pois- 
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Figure 9. Mass fraction in compact dark halo objects as a function of the mass of the lensing objects for combined LMC and SMC 
results from OGLE— II. The red box with a star shows the hypothetical value for the two LMC and one SMC events if they were caused 
by halo dark lenses. Solid lines show upper limits assuming all OGLE-II events are not due to MACHOs and the expected self-lcnsing 
background yields 3 events. Black dot-dashed lines show the same upper limit, but excluding the ambiguous SMC event. Also shown are 
the results of the MACHO collaboration (dotted line) and the upper limit derived by the EROS group (dot-dashed blue line). 



son statistics tables for such a background from Feldman & 
Cousins (1998). For a "typical" MACHO mass of 0.4 M© we 
expect at least 19 events, which gives a limit of 27 per cent 
of the total halo mass at 95 per cent confidence. 

However, if, based on its ambiguity, we reject the 
OGLE-SMC-01 as not a genuine microlensing event, there 
are only 2 events in both Clouds. The exact value of the 
background (i.e., self-lensing and other-than-MACHO lens- 
ing events) is not well known towards the SMC. Because, 
however, events detected by MACHO in the 1997 and 1998 
were confirmed in having a self-lensing nature, here we as- 
sume there should be at least one self-lensing event expected. 
For B = 3 and N = 2 the upper limit on MACHO fraction 
in the halo becomes about 20 per cent for deflector mass of 
0.4 Mq and 11 per cent for masses in range 0.003 and 0.2 
Mq , at the 95 per cent confidence level. The upper limit on 
the MACHO mass fraction is shown 4 in Fig. 9. 



9 CONCLUSIONS 

OGLE-II has provided a new and independent constraint on 
the presence of compact dark matter objects in the Galactic 
halo. In the LMC there were 2 candidate events found, both 
most likely due to self-lensing. A single candidate event was 
detected in the SMC data and its presence, if of microlens- 
ing nature at all, is consistent with the self-lensing scenario, 
in which the source is located at the far back end of the 
SMC and was lensed by a lens from within the SMC. The 

4 Data in electronic form are available through OGLE's website: 
http:/ / ogle, astrouw. edu.pl 



unusual position of the source on the colour-magnitude dia- 
gram could also be explained by, e.g., binarity of the source, 
and is generally not very trustworthy due to numerous am- 
biguities of the data. 

The derived optical depth estimate for the single event 
indicates a value of tsmc = (1.55 ± 1.55) x 10 -7 , which is 
very close to the previous measurements obtained with other 
data sets by the EROS and MACHO collaborations and is 
in agreement with self-lensing estimates. However, its low 
statistical significance prevents deriving any reasonable con- 
clusions on its origin, which may require further and more 
detailed studies of the structure of the SMC and resulting 
self-lensing optical depth. 

The detection of a single and unconvincing candidate 
event in the OGLE-II SMC data only strengthen our previ- 
ous conclusions reached in the OGLE-II LMC data (Paper 
I). The hypothesis that the Galactic halo is composed of 
compact objects is not favoured by these results. They show 
that the fraction of compact objects is close to zero and the 
only reason why we cannot completely exclude the compact- 
object hypothesis is due to limitations of the survey. 

This verdict, as determined from the OGLE-II data, 
will be further validated by analysing data from the recently- 
completed OGLE-III survey, which covers a much wider 
area and has a duration of 8 years. 
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